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Following HERA experiments our interest is focused, in the present 
paper, on charmonium production in lepton-proton collisions. 
Inelastic J/ ip electroproduction is quite difficult to be measured due 
to the low rates of events; almost all experimental studies are indeed 
performed in photoproduction regime. 

J/ ip resonance is considered as the product of the hadronization of 
a cc pair generated via boson-gluon fusion. 

The leading order contribution to the differential cross section, with 
respect to virtual photon-proton center of mass energy W and trans- 
ferred four-momentum squared Q^, is written as the sum of two 
terms related to a colour singlet and a colour octet cc pair produc- 
tion respectively. 

This expression contains two-gluons contributions and depends on 
three parameters; two of them are related to the hadronization frac- 
tions in charmonium states for the cc pair (colour singlet and colour 
octet), while the third parameter takes into accout possible "hot- 
spot" effects. 

The values of these parameters are estimated by mean of a fit to the 
experimental data collected by HI collaboration at HERA. 
By using different gluon densities, different sets of parameters are 
obtained, always with good fits. 
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Hidden charm production in e-p collisions at HERA is an intrigu- 
ing subject both from the theoretical and from experimental point 
of view. 

Available experimental data concern mainly charmonium photopro- 
duction which corresponds to limit (g^ = —Q^ is the trans- 

ferred four- momentum squared) [Q. 

In view of the low rates of events (cross sections decreasing with 
respect to Q^) only one sample of data (taken by HI detector) on 
J/ip electroproduction is available 0. 

It should be pointed out that e p e J/ip X photoproduction cross 
sections are characterized by one perturbative mass scale, namely 
the charm quark mass rric. 

These cross sections are found to be dominated by diffractive pro- 
cesses whose hallmark is the low value for X mass Mx- 
Coming to electroproduction events there are two perturbative mass 
scales which come into play, rric and Q. 

Either diffractive photoproduction or electroproduction experimen- 
tal rates have been recently understood in the framework of pertur- 
bative QCD (pQCD) i §. 



The main aim of this paper is the analysis of inelastic J/ip elec- 
troproduction; the true inelastic events are characterized by a Mx 
value much higher than the diffractive ones. 

The inelastic contributions to the inclusive cross sections can be 
therefore isolated by mean of a cutoff M^m on the values for Mx 

(HI coll. suggests Mmin = 10^ §). 

HI experimental points cover the kinematical range 

2 < g2 < 80 GeV^ 

so that electroweak corrections {Zq exchange) can be neglected and 
one can work in single-photon-exchange approximation. 
HERA results on open charm production in e~^p DIS suggest to 
describe also inelastic charmonium electroproduction in pQCD by 
cc pairs formation via boson-gluon fusion (BGF) mechanism given 
by the following partonic subprocess : 

■J* g ^ c c X 

That is to say any intrinsic charm contribution to the proton wave- 
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function is neglected. 




7* 

Figure 1: BGF mechanism 
The Born contribution is given by the following partonic diagrams: 




Figure 2: Born terms for e g ^ e c c 



Of course these diagrams contribute to the cross sections for pro- 
duction of a cc pair in a colour octet state (this pair is coupled to a 
gluon-photon system) 

e p — >• e |cc)g X 

The emission of an additional gluon allows the cc pair to emerge in 
a colour singlet state; so the lowest order contribution to the 

e p ^ e \cc)^ X 
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process comes from the following diagrams 




Figure 3: Lowest order contributions to 

Diagrams in Figure 3 contribute both to e ; 
e p — i> e |cc)g X. 



It should be pointed out that the two other diagrams shown in figure 
4 contribute to e p — > e |cc)g X process at the same perturbative 
order. 




Figure 4: Lowest order contributions to e g ^ e g |cc)g 



If we examine the overall lowest order contribution to the squared 
matrix element for e g ^ e g c c transition, we easily find that it 
enjoys the following colour decomposition: 

J2 \Mlo {eg^egcc)f = SqedT (F) C2 {F) + SnaC2 {A) (1) 

where the symbol ^ means that we are averaging over colour and 
polarization indexes for incoming particles and summing over out- 
going states. 

Sqed and Sna are the "space-time" coefficients; Sqed takes contri- 
bution only from the diagrams in figure 3, while S'tva includes also 
diagrams in figure 4. 

The colour factors are given, as usual, by 

_ 1 4 1 

C2(F) = ^^ = -, T(F) = -, C,{A) = N = Z 

The lowest order terms of the gluon initiated cross sections for a 
colour singlet cc pair inclusive production (e p — e |cc)^ X) turn 
out to be collinear finite since the charm quark is massive. 
Moreover the soft divergences found in the perturbative expansions 
for ag {e p ^ e |cc)g X) and ag {e p ^ e \cc)-^ X) compensate each 
other (order by order in perturbation theory); then the cross sections 
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(jg (e p — > e c c X) summed over colour states for the cc pair remain 
finite. They have been exphcitly computed in the hterature, up to 
next to leading order [^] |]IU| JTI] . 



If we consider the colour structure of the lowest order squared matrix 
element for the transition 

e g ^ e \cc)^ g 

we find it's given by 



lo (e 9 - e |cc), = SqedCi , Ci = (F) (F) 



where Ci is the corresponding colour factor. 

It's important to notice that HI events run over the ranges 

AO<W < 180 GeV , 0.2 < 2 < 1 

where W is the center of mass energy of the virtual photon-proton 
system and z is the ratio of J/ip and exhanged photon energies in 
the proton rest frame (elasticity of the process). 
It can be easily shown that the Bjorken scaling variable x, in terms 
of W and Q^, is given by x = ^ttj^- 

Accordingly the experimental points lie in the low x region 

6 ■ 10~^ < X < 5 • 10^2 

It's well known that the density of gluons pg (^,yu), with respect to 
the fraction of longitudinal momentum carried (^) and the factor- 
ization scale (fi), shows a rapid enhancement as ^ decreases. 
For X values small enough gluon recombination effects can no longer 
be neglected; in other words the BGF dynamics could involve 2 glu- 
ons (boson-gluon-gluon fusion). 

It should be stressed that these effects contribute to the cross sec- 
tions for inclusive colour singlet |cc)^ pair production already at 
the next-to-leading order; so we have explicitly calculated (under 
appropriate hypotheses) the lowest order contribution of the 

e g g ^ e \cc)^ 

partonic subprocess to the double-differential cross section with re- 
spect to W and variables 

(Per {e p ^ e \cc)^ X) 
dWdQ^ 
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The three partonic diagrams shown in figure 5 should be taken into 
account 




Figure 5: Boson-gluon-gluon fusion 

and those obtained by reversing the orientation of the charm 
quark fine. 6 diagrams in total. 

Their contribution can be written as a convolution over the two- 
gluons density for unit area in the proton with respect to the frac- 

(2) 

tions of longitudinal momentum carried Pg (^i, ^2, p): 

^2^(^o)(ep^e \cc), X) 
dWdQ^ 

j dCi j d^2 pf (ei, 6, P) F2, [6, 6, Q', S, W, as (p)] (2) 
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where y/S is the available center of mass energy for the electron- 
proton collision (at HERA ~ 300 GeV) and as (/u) is the strong 
running coupling constant. 

In the low density limit pg (^1,^2, A*) can be written as a product 
of single-gluon densities (the branches of a partonic cascade evolve 
independently) : 

Pf^ (Ci, 6, /^) = ^^P9 (^1' P) Pg (6, P) 

where R is the radius characteristic of the region in the proton po- 
pulated by gluons. 

It's natural to identify R with the proton radius R ~ Ifm (ghions 
are uniformely distributed inside the proton), but in the so-called 
"hot spot" effects R can be smaller (e.g R could be the radius of a 
constituent quark R < 0.5/m). 

It should be recalled that eq. (2) is infrared divergent because of 
soft singularities that occur in the limit Ci or ^2 — ^ 0. 
However gluon recombination favours gluons that are close in the ra- 
pidity space (they have greater probability to overlap and interact) . 
We make therefore the following change of variables 

[f] is the difference in rapidity for the gluon pair). 
Then eq. (2) becomes 

Pg (e-^e, p)CM) F,g [e, V, Q\ S, W, as (p)] (3) 

C (|?7|) is here a correlation function which implements mathemati- 
cally the previous guess, while 



Let us assume now that 

C{\v\)^S{rj) 
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This is, of course, a very strong hypothesis; it means that we allow 
recombination only for gluons with the same rapidity. Then one 
gets 

(Pa^i°\ep^e |cc), X) _ 
dWdQ^ 



/ (e, /^) [e, V = 0, Q\ S, W, as (/.)] (4) 

(0) 

Function ''7 = 0, Q^, S, W, as (a*)] in eq. (4) has been computed 

in pQCD, neglecting the transverse components of incoming gluons 
momenta and their interactions (free gluon wave functions). 
Its explicit expression is rather awkward, so we decided to put it in 
the appendix. 

Eq. (4) has been numerically integrated by using gluon densities 
implemented in CERN PDFLIB hbrary. 

The result turns out to be infrared safe; this is due to the hypothesis 
made on C (|7y|) shape which suppresses the initial state soft diver- 
gences. 

Of course this contribution has to be summed to the corresponding 
one-gluon term 

Therefore the leading order contribution to the hadronic differential 
cross section for e p — > e J/ip{z > zq) X process can be written as 
follows 

dV(^°) {ep^e J/jjj {z > zo) X) 



dWdQ^ 
d^(j{LO) (e p ^ e |cc)i X 



dWdQ'^ 

dV(^^) {ep^e \cc)^ X) 



n{\cc)^ J/iPiz > Zo)) + 
■n{\cc),^J/^{z>Zo)) (5) 



dWdQ^ 
with 

^2^(LO)(g^^g |cc)i X) _ 
dWdQ^ ~ 

dVf^^(ep^e \cc), X) d''a'if\ep^e \cc)^ X) 
dWdQ'^ ^ dWdQ'^ 
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The term 

£^iLO) (e p _ e |cc)8 X) 
dWdQ^ 

is the Born term which contributes to the production of a colour 
octet cc pair and 7i ^|cc)^ g — > J/ip {z > zo)j stand for the fraction 
of partonic events of type 

e p — >• e |cc)^ 8 X 

which hadronize in a J/ip vector meson with elasticity z > zq. 

As we have already mentioned — - — \mr^Qi — - — - is an infrared 
divergent quantity (soft singularities); on this subject, the next-to- 
leading order interference terms of the Born contribution with the 
virtual corrections for the e g ^ e c c transition (which contribute 
to a colour octet cc pair production) should be considered. 
In figure 6 arc shown some of them as cut diagrams. 
If we examine the colour structure of the interference terms we find 
that their contribution to the squared matrix element can be written 
as: 

J]/ (e ^ ^ e c c) = IqedT (F) (F) + Ij,aC2 {A) 

where Iqed, as before, takes contribution from "electrodynamical" 
diagrams (only quark-quark-gluon vertices), while In a contains the 
contributions of "non abelian" diagrams. 

Notice that (e g^e cc) has the same colour decomposition as 

^ \Mlo 9 ^ g 9 c both contributions to ^^-^ dWdQ'^ ^ c x) 

exhibit soft divergences; nevertheless the last quantity remains "in- 
frared safe" due to a compensation mechanism. 

This fact suggests to redefine — ^dWdQ'^^'^'^'^^ including part 

of the interference terms, which enter in — owdQ^ ' — ^' 

"^"'^ ^dWdQ^ c c X) (jgfing(j i-,y substituting 

\Mlo {e g^ e \cc)^ g)]"^ = SqedCi 

with 

Im}^'^^ {e g g cc) = {Sqed + Iqed) Ci 
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m the expression tor dWdQ^ ■ 

This operation has to be readsorbed in a redefinition of the hadroni- 
zation fractions Ti [\cc)^ g — > J /ip {z > zq)^. 









Figure 6: e g ^ e c c interference terms 



d'^a'-^^'>(e p^e c c X) . 



dWdQ'^ 



is found to be infrared safe in view of a Bloch- 
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Nordsieck compensation analogously to what happens in quantum 
electrodynamics; furthermore it can be obtained from the next-to- 
leading order contribution to '^^''^dwdo^ (computed by Smith 



et al. P |jTO[ [|TT|]) simply by fixing C2 (A) = and multiplying by 



the colour factor |. 
Then we add the two-gluons contribution and we find that 

^2~iLO) (^g p ^ g c c X) 

dWdQ^ ~ 

(P'd^^i ^ {e p ^ e c c X) (P'<y^2g'^^ {e p e \cc)-^ X) 
dWdQ^ ^ dWdQ^ 

Now eq. (5) can be written as follows 

dWdQ^ ~ 

^2~iLO) (^g p ^ g g g X) 



dWdQ"^ 
(i^af °) (ep^eccX) 



■Hi (cc-> J/^{z > zo)) + 
2 -Hsicc^ J/ij{z> zo)) ; (6) 



dWdQ 
where 

rf'gf (ep^eccX) _ d^a^^^^ {e p e |cc)g X) 
dWdQ^ ~ dWdQ^ 



We have already pointed out that the soft divergences in the per- 

' |cc>i x) 

compeni 

d?a[e p— >e |cc)g 



turbative expansion for *^ '^^'^ ^dWdQ^"^^ ^ ' compensate, order by order 



with analogous singularities in dWdQ^ ■ 

So the procedure introduced above, consisting in a subtraction of 

d^a[e p—i-e |cc)o X) d'^crle p— >e |cc), X) , 

terms from ^^^^^ ^ to ^vydQ^ ^' o^^er to com- 
pensate their soft divergences, can be extended, in principle, to all 
orders. 

In this way are defined the higher order ( "infrared safe" ) corrections 

, d^eri(e p^e c c X) i d'^ag,{e p^e c c X) 
dWdQ'2 '^'■^^ dWdQ^ 

Both these quantities depend, in general, on the renormalization 
scale n (which we set equal to the factorization scale), while the 
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physical cross section "'(^ '^~^dWdQ^^^°' depend on it. 

As a consequence 7ii and Hs have to depend in turns on /i in addi- 
tion to zq, W, and rric. 

In the following only zq dependence will be considered and all other 
effects will be neglected. 

HI collaboration has experimentally estimated the values of the dif- 

r , . 1 , . da(e p^e J/ip(z>0.2) X) i da(e p^e J/i>(z>0.2) X) 

ferential cross sections — ^-^ — - and — ^-^ — - 

for inelastic J/ip electroproduction in p collisions. 
Therefore we have numerically integrated our leading order expres- 
sion in eq. (6), with respect to Q'^ or W. 
The factorization/renormahzation scale is set as follows 

The right hand side of eq. (6) is then fitted to experimental data by 
using the minirnization routines in MINUIT package of F.James 
and living Hi, lis and R as free parameters. 

Different choices for the density of gluons are made, in order to com- 
pare the fit's results. 

The following NLO gluon densities have been used: MRS (G) 
(02.1995), CTEQ 2pM, CTEQ 4M and GRV 94 HO. 
The results of these calculations are reported in table 1. 
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MRS (G) (02.1995) gluon density 




R = 


-- l.lO/m , Hi (^0 = 0.2) = 0.269 , Tig (-2o = 0.2) = 
= 6.74 , d.o.f = 10 ^ Xred = 0.674 


= 2.75- 


10~3 




CTEQ 2pM gluon density 






i? = 


0.442/m , Til (zo = 0.2) = 0.175 ,7Yg(^o = 0.2) 
X' = 6.32 , d.o.f = 10 ^ xfed = 0.632 


= 4.21 


•10-3 




CTEQ 4M gluon density 








0.393/m , Hi {zq = 0.2) = 0.109 , Tig {zq = 0.2) 
X" = 6.23 , d.o.f = 10 ^ = 0.623 


= 4.79-10-3 




GRV HO gluon density 






R = 


0.611/m , ni{zo = 0.2) = 0.160 , Tig (;So = 0.2) 
X' = 6.51 , rf.o./ = 10 ^ xled = 0.651 


= 5.32 


•10-3 



Table 1: Comparison of fit results for different gluon densities 
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In the following graphs HI experimental points are compared 
with our results. 

MRS (G) (02.1995) gluon density 
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CTEQ 2pM gluon density 
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CTEQ 4M gluon density 
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GRV HO gluon density 
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The errorbars are built-up summing in quadrature statistical and 
systematical errors (obviously has been computed taking into ac- 
count only statistical ones). 

We remark that, even if these fits give values for Tig 1-2 orders of 
magnitude smaller than for TYi, the perturbative terms connected to 
Til are one order of magnitude smaller than those linked with Tig; 
the latter contribution can be up to one order of magnitude smaller 
than that associated with Hi . 

The large discrepancies in the fitted parameters obtained using dif- 
ferent gluon densities indicate that inelastic charmonium production 
provides a particularly sensitive tool to probe the gluonic structure 
of the hadrons; this is due to the fact that the two-gluons contri- 
bution, which depends on the square of the gluon density and then 
varies strongly with the choice of the density itself, is of the same 
order of magnitude as the single-gluon contribution (of the same 
perturbative order) in our kinematical range. 

This fact suggests finally that important informations on the dis- 
tribution of gluons in the transverse space in the proton could be 
obtained from hidden charm electroproduction. 
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Appendix 

1 



1 2 



Ami 



{ 



1920m^ 



16m' 



(2^-1) + 21^2^ 

2^2 _ 25 (g2 + w^) + (g2 + w^y] + 

86 + V^2^) _ 2^ ^g2 ^ ^,2^) ^ (g2 ^ ^,2^)2j ^| ^ 

2>S' - 25 + 1^2^) ^ (^g2 ^ ^2^)21 _ 1) + 2W^(\ 



+ 2m; 



{-32Q^ 



+ 9605Q^ 

2 .2 



-165^ + 165 (g^ + W^) + (Q^ + W^f 

(q2 + (g^ + - 5) e + 44 (g2 + vr^) ' e 

252 - 25 + W'^) + {Q'^ + W'^Y } } } + {320m^ 
252 - 25 + w'') + + w^2^2l ^ |^g^2 

552 - 55 + + + ly^)'] - 42 (g2 + 

252 - 25 + V^2^ ^ ^q2 ^ ^2^2j ^1 ^ 2^^2 (5 _ g2 _ p^2^ 



32Q^ - 112Q2 (Q2 + |y2^ ^ ^ ^g2 ^ 1^2^ 2 ^2 

{2565g^ (5 - - w^') - 56g=^ (g2 + w^) [105^- 
105 (g2 + w^) + (g2 + w^f] c + 84 (g2 + w'Y \{q' + w'Y + 



25^-25 {Q' + W')] e}}bog 



1 + 



4m2 



Q2(25-1)+2W^25 



1 - 



4m2 

Q2(2^_1)+2W2$ 



(7) 
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where a ~ is the e.m fine structure constant and Qc — ^ repre- 
sents the quark charm electric charge (in e units). 
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